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Recent work on the numerical solution of supcrsymmctric gauge theories is de- 
scribed. The method used is SDLCQ (supersymmetric discrete light-cone quan- 
tization). An application to N = 1 supersymmetric Yang— Mills theory in 2+1 
dimensions at large N c is summarized. The addition of a Chern-Simons term is 
also discussed. 



1. Introduction 

Although much has been learned about supersymmetric gauge theories by 
analytic methods, numerical methods can yield much more of the nonper- 
turbative structure. In particular, the method known as supersymmet- 
ric discrete light-cone quantization (SDLCQ), 1,2 an extension of ordinary 
DLCQ, 3,4 has been quite successful in the analysis of (l+l)-dimensional 
supersymmetric theories. This work has recently been extended to 2+1 
dimensions 5,6,7,8 with consideration of TV = 1 supersymmetric Yang-Mills 
(SYM) theory, including a Chern-Simons (CS) term. 9 The mass spectrum, 
Fock-state wave functions, and a stress-energy correlator are all computed. 
The CS term provides an effective mass that reduces the tendency of SYM 
to produce stringy states with many constituents. This work was done at 
large- N c , but the method is also applicable to finite N c . 

As the name SDLCQ implies, light-cone coordinates 10 are used. They 
are defined by spacetime coordinates 

x ± = (t±z)/V2, x x = (x,y) (1) 
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and momentum components 

p ± = (E± Pz )/V2, p ± = (p x , Py ). (2) 

The dot product of two such four- vectors then becomes 

p ■ x = p + x~ + p~x + — pj_ • xj_ . (3) 

The x + direction is treated as the direction of time evolution, which makes 
the conjugate variable p~ the light-cone energy. The light-cone three- 
momentum is p = (p + , pj_). In a frame where the net transverse momentum 
P_i_ is zero, the mass eigenvalue problem becomes 

2P+P-\P) = M 2 \P) , (4) 

where \P) is also an eigenstate of three-momentum P. The SDLCQ method 
provides a means to solve this eigenvalue problem with supersymmetry 
preserved exactly at any level of the approximation. 

One of the advantages of light-cone coordinates is that there exists a 
well-defined Fock-state expansion for each mass eigenstate. There are no 
disconnected vacuum pieces, because the longitudinal momentum of each 
constituent, virtual or real, must be positive. SDLCQ uses a Fock-state 
expansion for |P) to obtain a matrix eigenvalue problem for the Fock-state 
wave functions at discrete values of the momentum. The matrix is then 
diagonalized by appropriate means. For large matrices the Lanczos diago- 
nalization technique 11 has been used, as discussed in Ref. [ 6]. 

The discretization is accomplished by restricting the fields to periodic 
boundary conditions in a light-cone box 3 ' 4 defined by — L\\ < x~ < L\\ and 
< x, y < L±. This leads to a discrete momentum grid 

P ^>~rn, p± -» (— n x , — n y ) . (5) 

The product P + P~ is independent of Lp and the limit — > oo is ex- 
changed for a limit in terms of an integer K, called the harmonic resolution, 3 
defined by 

K = ^P+ . (6) 

Longitudinal momentum fractions x = p + /P + then reduce to n/K. The 
number of particles in a Fock state is limited to K, because negative longi- 
tudinal momentum is not allowed and the individual integers n must sum to 
K . a The eigenvalue equation (4) becomes a coupled set of integral equations 
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for the Fock-state wave functions in which the integrals are approximated 
by discrete sums over the momentum grid 

f dp+ I rfV/(p + ,Pi) * £~ (if) E f(nP + /K,2n ± ir/L ± ) . (7) 

The harmonic resolution provides a natural cutoff for n. The transverse 
sums must be truncated explicitly, which is done by limiting n x and n y to 
range from — T to T. The integer T can be viewed as a transverse cutoff or, 
at fixed dimensionful cutoff Aj_ = 2nT/L±, as the transverse resolution. 

The distinction between DLCQ and SDLCQ lies in the choice of operator 
for discretization. In ordinary DLCQ one discretizes the Hamiltonian, P~; 
in SDLCQ, one discretizes the supercharge Q~ and constructs P~ from the 
superalgebra relation 

{Q-,Q-} = 2V2P- , (8) 

which guarantees that the discrete eigenvalue problem preserves 
supersymmetry. 1,2 The P~ of ordinary DLCQ differs from the supersym- 
metric P~ by terms which disappear in the large- if limit but which break 
the supersymmetry at finite K . 

The remainder of this paper is organized as follows. In Sec. 2, (2+1)- 
dimcnsional SYM theory is reviewed and numerical results discussed for the 
spectrum and for a correlator of the stress-energy tensor. Section 3 extends 
the study of the spectrum to include the CS term, in both a dimcnsionally 
reduced theory and the full (2+l)-dimensional case. A brief summary is 
given in Sec. 4. 



2. SYM 2+ i theory 
2.1. Formulation 

The action for N = 1 SYM theory in 2+1 dimensions is 

S = J dx+dx-dx^ii—F^F^+i^D^), (9) 

with 

F liV = d li A v -d v A ll +ig[A ll ,A v ], D„ = + ig[A„ }. (10) 
The fermion field is separated into chiral projections 

1 + 1 - T 5 
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only one of which is dynamical. In light-cone gauge, A + = 0, with the 
transverse component of the gauge field A± written as <j>, the action becomes 



/ 



dx + dx dx±tr 



l -{d-A-f + D + t 



> + iipD + ip+ 



+ixd-x + ■j=i/>D±<j> + -^=(f>D_Ltp 



The non-dynamical fields A and \ satisfy constraint equations 

1 



(12) 



(13) 



by which they can be eliminated from the action. The dynamical fields are 
expanded in terms of creation operators 



| 00 PC 

V^Ll Jo 

n — — oo 

+a] i {k+,n ± y 



dk+ 
V2k+ 



aij(k + , n )e 



-ik + x — i- 



(14) 



ipij{0,x ,x±) = 



E 



dk + 



b i j(k + ,n ± )e 



-ik~*~ X — i- 



+bUk+,n^e ik+x - +il 



(15) 



where in (2+1) dimensions n 1 - is the only transverse momentum index. 
The chiral components of the supercharge are 



d^W + g (i[0, 8-4>] + 2^) 



(16) 
■ (17) 



Q+ yM = 2 1 / 4 J dx~dx ± tv [cj)d-ip - ipd-<j)} , 

Qsym ~ 2 3 ^ 4 J dx~dx±iv 
They satisfy the supersymmetry algebra 
{Q + ,Q + } = 2V2P+, {Q-,Q-} = 2V2P- , {Q + ,Q-} = -4P± . (18) 

This theory has the additional symmetries of transverse parity, P: 
aij(k,n ± ) — ► —a i j(k,—n ± ), b i j(k,n ± ) — > b i j(k,—n ± ) and Kutasov's 
transposition 12 S: aij(k,n J -) — » —aji(k,n- L ), bij(k,n J -) — > —bji(k,n ± ). 
These allow the matrix representation to be block diagonalized by an ap- 
propriate choice of basis. Eigenstates are labeled by the quantum numbers 
±1 associated with P and 5. 
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2.2. Spectrum and wave functions 

The main results for the spectrum of the SYM 2 +i theory are given in Figs. 1, 
11, and 12 of Ref. [ 6]. They show that the masses squared can be classified 
according to three main forms of behavior: 1/L\, g 2 N c A±, and A\. In 
particular, the spectrum as a function of g separates into two bands, one 
of approximately constant M 2 L\ and the other growing rapidly. 

For states in the lower band, the average number of constituents in- 
creases rapidly with g. At g ~ 1.5W4n 3 /N c L± the DLCQ limit of K 
constituents is saturated. Thus in SYM theory the low-mass states are 
dominated by Fock states with many constituents, in close correspondence 
with string theory. However, as a practical matter, the saturation means 
that the SDLCQ approximation breaks down, and numerical studies in this 
band must be limited to smaller couplings. 

Within the coupling limitation, extrapolations to infinite resolution are 
easily done for low-mass states. One first considers M 2 as a function of 1/T 
for a sequence of fixed K values. The extrapolations to T = oo then yield 
M 2 as a function oil/K alone, to extrapolate to K — oo. The different 
representatives of continuum eigenstates are disentangled by studying their 
properties, such as average constituent content and momentum. Of course, 
the different P and S symmetry sectors are explicitly separated at the start. 
Typical extrapolations are illustrated in Ref. [ 6] with plots in Figs. 4-8 and 
results in Tables II and III. 

For the spectrum as a whole there is a curious behavior with respect to 
the average number of fermion constituents (rip). Calculations for trans- 
verse resolution T = 1 and longitudinal resolutions K = 5 and 6, and for 
many different coupling strengths, show a gap between (rip) =4 and 6, 
where no state is found. 

Wave functions are also obtained in the diagonalization process. In 
the analysis of the spectrum they were used to compute various average 
quantities that helped identify states computed at different resolutions. 
More of the form of the wave function is revealed in the structure function 

K K-q T / q \ / q 

g A (n,n^) = J2 £ £ sl^m-KU £n+ 

qr=2 ni,...,n, = l n^,...,n^-=-T \i=l / = 1 

x £ W8$ 6 i \^nh • • • ; n q , <)| 2 , (19) 
i=i 

where A and A\ represent the statistics (bosonic or fcrmionic) of the probed 
type and the i-th constituent, respectively, and tp is a Fock-state wave func- 
tion. In the lower band the shapes are typically simple and are found to 
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confirm the identification of states at different resolutions. In the upper 
band, there are complicated shapes with multiple bumps in transverse mo- 
mentum, such as in Figs. 13 and 14 of Ref. [ 6]. 



2.3. Stress- energy correlator 

Consider the following correlator of the stress-energy component T ++ : 

F(x+,x-,x ± ) = (0|T ++ (x+,a;-,a; i )r ++ (0,0,0)|0) , (20) 

at strong coupling, 7 as an example of what one might compare with a 
supergravity approximation to string theory for small curvature. 13 In the 
discrete approximation, F can be written as 

F(x+,x-,0)= £ (") (21) 



where 



n,m,s,t 

x {0\^-T(n, m)e- lP °? x+ - lP+x ~ ^-T(s, t)\0) , 

7T 7T 



^T++(n,m)|0) = ^tr [at.(n, nx)at(m, m ± j\ |0) (22) 



(n — m) 

- tr 



bl j (n,n±)b] i (m,m±) |0) 



4 

Insertion of a complete set of bound states \a) with light-cone energies 
P~ = {Ml + Pj)/P + at resolution K (and therefore P+ = itK/L\\) and 
with total transverse momentum P± = 2Ttt/L^ yields 

1 f X ~V „t + ^^ V- 1 M l' 2 r, „„ \ \(u\a)\ 2 

(23) 

x± = 0, and 

|u) = N u ^-y2d n+m , K S n±+m±tN± T(n,m)\0) . (24) 

IT z — ' 

n.m 

Here iV„ is a normalization factor such that (u\u) = 1. The sum over the 
full set of eigenvalues can be avoided by a Lanczos-based technique. 14 

For free particles, (x + /x~) 2 F has a 1/r 6 behavior. 7 In the interacting 
case, this behavior should be recovered for small r, where the bound states 
behave as free particles. Because this behavior depends on having a mass 
spectrum that extends to infinity, the finite resolution of the numerical 



7 



calculation yields only 1/r 5 ; however, the 1/r 6 behavior is recovered in a 
careful limiting process. 

The behavior for large r is determined by the massless states. Because 
this theory has zero central charge, there are exactly massless Bogomol'nyi- 
Prasad-Sommerficld (BPS) states at any coupling. However, their wave 
functions remain sensitive to the coupling, and, at a particular (resolution 
dependent) value of g, the correlator is exactly zero in the large-r limit. 
The associated 'critical' value of g increases in proportion to the square 
root of the transverse resolution T. 7 



3. SYM-CS theory 

The following CS term can be added to the Lagrangian: 

L CS = \^ vX (a^Ax + j 9 A„A v A x ^ + . (25) 

This induces an additional term in the supercharge 

kQcs = -2 3/4 «y dx-d-cf)-^-iP , (26) 

which generates in P~ terms proportional to k 2 that act like a constituent 
mass squared. The presence of an effective mass reduces the tendency for 
low-mass states to be composed of high-multiplicity Fock states. This cre- 
ates a theory in which the eigenstates are more likely to be QCD-like, i.e. 
valence dominated, and improves the applicability of the SDLCQ approxi- 
mation to a greater range of couplings. 15 ' 16 ' 8 

The dominance of the valence state is most prominent in the dimcn- 
sionally reduced theory. 15 This (l+l)-dimensional theory is obtained by 
requiring the fields to be constant in the transverse direction and replacing 
d± by zero in the full supercharge <5g Y M + K Qcs- The SYM contribution 
is then proportional to g. Figure 3 of Ref. [ 15] illustrates the dramatic re- 
duction in the average number of constituents as the ratio n/g is increased. 
Also, there are states for which the mass is nearly independent of g at fixed 
k. These are identified as approximate BPS states and are the reflection 
of the massless BPS states in the underlying SYM theory. 16 This behavior 
can be seen in Fig. 1 of Ref. [16]. 

Similar anomalously light states appear in the full (2+l)-dimensional 
theory. 8 The bulk of the spectrum is driven to large M 2 values as g is in- 
creased, but one or more states remain at low values. The presence of the 
transverse degree of freedom makes this more difficult to disentangle numer- 
ically, because the matrices are larger and because one must consider the 
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transverse resolution limit. Also, the eigenstates are less valence-dominated 
at stronger YM coupling, which makes the SDLCQ approximation less use- 
ful. However, structure functions have been extracted at intermediate cou- 
pling to show that the approximate BPS states have a distinctly flat depen- 
dence in longitudinal momentum. Figure 3b of Ref. [ 8] gives an example 
of this behavior. 

4. Summary 

This work shows that one can compute spectra, wave functions, and matrix 
elements nonperturbatively in supersymmctric theories. The introduction 
of a Chern-Simons term brings an effective constituent mass which has the 
effect of reducing the tendency of SYM theory to form stringy, low-mass 
states with many constituents. Instead, the lowest-mass states tend to be 
dominated by their valence Fock state. The massless BPS states of SYM 
theory survive in SYM-CS theory as states with masses nearly independent 
of the YM coupling. 

A number of interesting issues remain to be explored. They include 
theories in 3+1 dimensions, matter in the fundamental representation, 17 
and supersymmetry breaking. All of these are important for making contact 
with QCD. 
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